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I.  INTRODUCTION 

A.  THE  MARGINAL  ICE  ZONE  AND  THE  EAST  GREENLAND  CURRENT 

The  marginal  ice  zone  is  a  region  where  complex  atmospheric  and  oceanic 
physical  processes  affect  the  distribution  of  ice  between  the  pack  ice  and  the  open 
ocean.  In  recent  years  the  marginal  ice  zone  (MIZ)  in  the  region  of  the  East 
Greenland  Current  has  been  extensively  investigated  during  Marginal  Ice  Zone 
Experiments.  The  East  Greenland  Current  is  arctic  (polar)  water  flowing  southward 
along  the  east  coast  of  Greenland,  over  the  200  m  deep  and  200  km  wide  continental 
shelf,  and  then  continuing  to  flow  southward  over  deeper  (4000  m)  water  (Wadhams  et 
ai,  1979).  The  boundary  formed  by  the  arctic  water  of  the  East  Greenland  Current 
and  the  Atlantic  water  is  a  potentially  unstable  ocean  front.  Figure  1.1  shows  the 
location  of  the  East  Greenland  Current  and  the  mean  ice  edge. 

The  study  of  stability  of  ocean  fronts  has  been  done  in  the  laboratory  as  well  as 
in  numerical  experiments.  Griffiths  and  Linden  (1981)  examined  stability 
characteristics  of  density  driven  fronts  in  rotating  tank  experiments.  Their  experiments 
indicate  that  the  type  of  instability  (barotropic  or  baroclinic)  associated  with  a  given 
flow  field  is  dependent  upon  the  length  scale.  In  general,  if  the  length  scale  L  of  the 
mean  flow  is  larger  than  the  local  internal  Rossby  radius  of  deformation  (Rj),  then  the 
resulting  instability  is  baroclinic.  Where  the  length  scale  is  comparable  to  R^,  the 
instability  is  more  likely  to  be  barotropic.  The  width  scale  of  the  East  Greenland 
Current  is  on  the  order  of  200  km,  much  greater  than  the  local  Rj  (8  km),  suggesting 
baroclinic  instablity  (Griffiths  and  Linden,  1981).  The  observed  50  km  length  scales  of 
some  of  the  eddies  (Wadhams  and  Squire,  1983)  appear  to  be  consistent  with  these 
considerations.  The  mesoscale  ice  edge  frontal  eddies  observed  by  Johannessen  et  al. 
(1983)  had  a  length  scale  comparable  to  the  local  internal  Rossby  radius;  therefore, 
from  the  above  considerations,  these  mesoscale  eddies  would  be  the  result  of  barotropic 
instability. 

B.  EDDIES  IN  THE  EAST  GREENLAND  CURRENT  REGION 

From  an  observational  standpoint,  eddies  in  the  MIZ  have  only  recently  been 
documented.  These  eddies  have  been  observed  to  have  scales  ranging  from  5  to  80  km. 
There  are  five  possible  sources  for  the  eddies  of  the  East  Greenland  Current  region: 
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Figure  1.1     Ice-Edge  Oceanic  Front  (from  Johannessen  et  al,  1987). 


(1)  Barotropic  instability  was  suggested  by  Johannessen  et  al.  (1983)  as  being 
responsible  for  the  formation  of  small-scale  eddies;  (2)  Baroclinic  instability  is  largely 
responsible  for  the  larger-scale  eddies  (Griffiths  and  Linden,  1981);  (3)  Current- 
topography  interaction  may  cause  eddy  formation  (Smith  et  al,  1984);  (4)  Open  ocean 
eddies  present  in  Atlantic  water  may  be  advected  toward  the  ice  edge,  leading  to 
interaction  that  can  develop  into  ice-edge  eddies  (Johannessen  et  al.,  1987);  and  (5)  Ice 
edge  air-sea  interaction  may  cause  eddies  to  form  (Reed  and  O'Brien,  1983; 
Hakkinen,  1986a,  b). 


C.  MESOSCALE  EDDIES 

Studies  of  small-scale  (5  to  15  km)  eddies  observed  in  conjunction  with  the  ice 
edge  front  over  deep  water  (Johannessen  et  ai,  1983  and  Johannessen  et  ai,  1987)  have 
significantly  increased  the  knowledge  of  the  structure  of  mesoscale  eddies.  The  front 
was  observed  to  have  velocities  of  approximately  10  cms  (with  little  vertical  shear)  and 
a  width  scale  of  10  km.  Observations  found  several  mesoscale  eddies  with  length  scales 
ranging  from  5  to  15  km  associated  with  the  front.  Like  the  front,  these  associated 
eddies  also  exhibited  little  vertical  shear.  Based  on  considerations  from  theory  and 
laboratory-  experiments  discussed  above,  Johannessen  et  ai  (1983)  reasoned  that  the 
eddies  most  likely  resulted  from  barotropic  instability  of  the  oceanic  front. 
Johannessen  et  ai  (1987)  also  established  that  eddies  and  meanders  are  the  dominant 
features  along  the  ice  edge  under  moderate  wind  conditions.  Mesoscale  eddies  are  the 
focus  of  this  study. 

D.  WIND-FORCED  STUDIES 

Numerous  studies  of  wind  driven  ocean  motion  in  the  marginal  ice  zone  have 
been  conducted  since  1980.  Initially  these  were  two-dimensional  cross-ice  edge 
examinations  of  upwelling  and  downwelling  produced  by  along-ice  edge  winds  (Roed 
and  O'Brien,  1983;  Hakkinen,  1986a).  Because  the  coupling  is  stronger  under  ice, 
Ekman  transport  is  larger  under  ice.  Roed  and  O'Brien  (1983)  describe  the  relationship 
between  the  wind  and  the  ice  edge  as  "...  winds  which  blow  along  the  ice  edge  and  to 
the  left  when  facing  the  ice  which  favors  upwelling."  Hakkinen  (1986b)  also  showed 
that  along  ice  edge.  10  m  s  winds  produce  upwelling  (downwelling)  when  the  ice  edge  is 
to  the  right  (left)  of  wind  direction. 

Recently,  Smith  et  ai  (1987)  have  studied  the  effect  of  along-ice  edge  winds  on 
existing  ocean  eddies.  Their  results  show  that  for  constant  10  m  s  winds,  the  ice 
responds  largely  to  the  wind  and,  to  a  lesser  extent,  the  ocean  eddy.  They  also  find, 
however,  that  upper  ocean  anticyclones  (cyclones)  are  rapidly  destroyed  by  upwelling 
(downwelling)  favorable  winds.  None  of  these  studies  have  considered  the  effects  of 
lighter  winds  and  reversals  in  the  wind  direction. 

Simulations  using  constant  wind  directions  are  not  altogether  realistic.  Figure  1.2 
and  Figure  1.3  show  the  wind  pattern  observed  by  Johannessen  et  ai  (1987)  and 
Morison  et  ai  (19S7),  respectively,  during  MIZEX  84.  The  two  observations  show 
that  there  were  periods  when  the  wind  shifted  regularly,  and  also  there  were  extended 
periods  with  light  winds. 
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Figure  1.2    Time  Series  of  Wind  (from  Johannessen  et  al.  (1987)). 


15.2  cm  f1 
i 1 


1  '  '  '  '  '  '  '  '  '  '  '  '  '  '  I  I  I  I  I  I  I  I  I  I  '  I  I 

160     1G«     ,68     »72     t7G     180     164     188     192     196     200 

DAY  OF  1964 

Wind  and  ice  veloc.ty  vector,  at  the  Polar  Qmn,  The  wind,  are  6-hour  average,  of  measured  wind  at  10  m.  centered  at 
the  time  ihown. 


Figure  1.3     Time  Series  of  Wind  (from  Morison  et  al.  (1987)). 
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Johannessen  et  al.  (1987)  provide  a  schematic  of  the  interaction  of  an  open 
ocean  eddy  with  the  ice  edge  (Figure  1.4).  This  eddy  (Johannessen  et  al.  (1987),  eddy 
(El 3))  was  observed  for  about  eleven  days.  Figure  1.2  shows  that  for  the  first  three 
days,  steady  winds  along  the  ice  were  observed,  followed  by  periods  of  lighter  winds, 
varying  in  direction  for  periods  of  about  one  day.  The  last  days  of  observation  of  this 
eddy  were  dominated  by  strong  winds  from  the  south,  forcing  the  ice  away  from  the 
eddy. 

E.       PURPOSE  OF  THIS  STUDY 

The  objective  of  this  study  is  to  understand  the  interaction  of  mesoscale,  and 
largely  barotropic,  eddies  with  a  marginal  ice  zone  through  the  use  of  a  nonlinear  two- 
layer  regional  ocean  model  with  temporally  varying  winds  and  a  coupled  ice  model. 
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Schematic  of  the  interaction  of  the  open  ocean  eddy  (E13) 
*uh  the  ice  edee 


Figure  1.4    Open  Ocean  Eddy  Interacting  with  Ice  Edge. 

From  Johannessen  et  al.  (1987).  Where  AW  is  Atlantic  Water,  MAW  is  MIZ  Atlantic 
Water,  ASW  is  Arctic  Slope  Water,  and  MASW  is  MIZ  Arctic  Slope  Water. 
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II.  NUMERICAL  TECHNIQUE  AND  MODEL  PARAMETERS 

A.       OCEAN  MODEL 

1.  Ocean  Model  Equations 

Experiments  are  performed  using  a  numerical  model.  The  model  is  based  on 
primitive  equations,  has  two  layers,  and  uses  a  semi-implicit  numerical  scheme. 
Motion  in  each  layer  is  governed  by  a  momentum  equation: 

dV-Jdt  +  (VV-  +  V-V)  v:  +  k    x  fv-  = 

1  l  l         l  l  (eqn2.1) 

-hlV?{  +  AhV2Vi  +  (8il/p1)((l-A)taw+tlw) 
and  a  continuity  equation: 

dhj/at  +  V  •  Vi  =  0  (eqn  2.2) 

for  layer  (i=  1  upper  and  i=  2  lower),  thickness  h-,  transports  V-,  and  velocities  vj.  The 
fluid  is  hydrostatic,  Boussinesq,  and  the  fluid  density  (p-)  in  each  immiscible  layer  is 
fixed.  Subgrid  scale  turbulent  eddy  dissipation  processes  are  represented  by  a 
horizontal  Laplacian  operator  on  transport.  As  this  only  approximates  a  complex 
subgrid  scale  dissipation  process,  its  coefficient,  A^,  has  been  chosen  small  (10  m  /s). 
The  fluid  is  thus  relatively  inviscid.  The  notation  used  in  all  equations  is  defined  in  the 
Appendix.  The  above  scheme  has  been  used  in  numerous  mid-latitude  ocean  mesoscale 
circulation  studies  (Hurlburt  and  Thompson  (1980,  1982);  Smith  and  O'Brien  (1983)) 
where  more  thorough  presentations  show  how  energy  and  enstrophy  are  conserved  in 
the  absence  of  dissipation.  Smith  and  Reid  (1982)  compared  test  cases  to  linear 
analytic  solutions  to  verify  the  Rossby  dispersion  characteristics  of  the  model. 

2.  Boundary  and  Initial  Conditions 

A  rectangular  (55  x  40  km)  finite  difference  gridded  domain  is  used.  There  are 
55  grid  points  in  the  x-direction  and  44  grid  points  in  the  y-direction.  The  domain  is 
rotated  counter-clockwise  approximately  30  degrees  relative  to  a  reference  latitude  of 
80  °N  so  that  the  x-axis  is  aligned  with  the  East  Greenland  Current,  approximately 
parallel  to  the  east  coast  of  Greenland.  The  initial  ocean  state  consists  of  a  Gaussian 
height  field  in  each  layer,  which  is  centered  in  the  basin: 
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h.  =  Aoe(V  +  y2)/2L2))  (eqn23) 

L  is  the  e-folding  width  scale  for  the  eddy  ( =  5  km).  The  depth  of  thermocline  in  the 
East  Greenland  Current  ranges  from  25  to  100  m.  All  experiments  for  this  study  are 
initialized  with  the  upper  layer  mean  thickness  of  50  m  and  the  lower  layer  mean 
thickness  of  4000  m.  The  resulting  first  internal  Rossby  radius  of  deformation  (R^)  is 
equal  to  5  km,  consistent  with  the  observations  of  Johannessen  et  al.  (1983). 
Nondimensional  eddy  size,  y=L/R(j,  is  equal  to  1.  The  amplitude  (AQ)  of  the 
Gaussian  distribution  was  chosen  so  that  a  maximum  velocity  would  be  on  the  order  of 
10  cm/s  in  each  layer.  Thus,  the  Rossby  number  (vmax/fL)  for  such  a  flow  is  0.15; 
therefore,  the  nonlinear  terms  are  expected  to  be  important  parameters  in  the 
momentum  equations.  The  ocean  velocity  field  is  thus  initially  in  gradient  balance  and 
constant  with  depth  (uj  =  U2).  The  north  and  south  boundaries  of  the  model  domain 
are  no-slip  walls,  where  both  tangential  flow  and  normal  flow  are  set  equal  to  zero. 
Open  radiation  is  allowed  at  the  east  and  west  boundaries  of  the  model  domain  to 
allow  wind-driven  inflow  and  outflow.  Simulations  were  integrated  for  short  periods, 
typically  4  days.  The  sense  of  ocean  eddy  rotation  (cyclonic  vs.  anticyclonic),  wind 
direction,  wind  speed,  and  the  duration  of  wind  shifts  were  varied.  Although  the  eddies 
are  initialized  with  barotropic  structure,  baroclinic  structure  associated  with  wind- 
driven  layer  interface  motion  is  free  to  evolve. 

B.       ICE  MODEL 

1.  Ice  Model  Equations 

Motion  in  the  ice  is  governed  by  momentum  equations: 

dujdx  +  uduldx  +  vdu/dy  =  fv  -(A/m)(Txai  -  Txiw)  -g*^(h1  +  h2)/5x  <ecln  2-4) 


dv/dt  +  udv/dx  +  vdv/dy  =-fu  -(A/m)(t  ai-  x  1W)  -%*d(hx  +  h2)/dy  (eqn  2.5) 

and  by  continuity  equations: 

3A/dt  +  d(Au)/ax  +  d(Av)  >dy  =  A„72A  ^n  2-6> 
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dm/dt  +  d(mu),'dx  +  d(mv)  fdy  =  AmV2m  (eqn  2.7) 

for  ice  concentration  (A)  and  ice  mass  (m).  As  in  previous  ice  model  studies,  a 
Laplacian  damping  term  has  been  included  in  the  continuity  equations  for  A  and  m. 
The  ice  equations  contain  a  pressure  force  associated  with  sea  surface  slope.  Most 
mesoscale  modeling  studies  of  the  MIZ  have  neglected  this  term;  Smith  et  al.  (1987) 
show  that  this  term  is  important  in  the  absence  of  wind  forcing.  No  internal  ice  stress 
was  used,  which  is  suitable  for  concentrations  less  than  85%.  Thermodynamic  effects 
are  not  considered  in  this  study.  The  ice  is  coupled  to  the  ocean  through  an  ice-water 
interfacial  stress,  T1W,  for  ice  (u)  and  ocean  (u^,)  velocity  vectors: 

T'W=  PiSw^-Mu-Uwl  <ecln  2-8> 

and  similarly  the  ice  is  coupled  to  the  air  through  the  air-ice  interfacial  stress,  Ta\  for 
air  (ua)  and  ice  (u-ce)  velocity  vectors: 

Tai=Pacai(Vuice)IVuicel  (eqn  2.9) 

All  constant  values  (p-,  pa,  c-w,  c^)  are  chosen  following  Hakkinen  (1986a)  and  are 
given  in  the  appendix. 

2.  Ice  Thickness  and  Concentration 

The  ice  thickness  distribution,  D,  is  initially  specified  to  be  2  meters  thick,  but 
is  then  allowed  to  vary  according  to 

D  =  m/(pjA)  (eqn  2.10) 

Ice  concentration  is  initially  a  concentrated  band  in  the  middle  of  the  domain.  The  ice 
band  is  a  linear  function  of  the  y-coordinate,  varying  from  0.05  at  the  edge  of  the  band 
nearest  the  open  ocean  to  0.75  at  the  edge  of  the  ice  band  nearest  the  pack  ice.  The 
ice  is  initially  at  rest.  The  difference  in  value  between  adjacent  contour  lines  of  ice 
concentration  found  in  the  model  output  graphics  is  0.05. 
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C.       PRELIMINARY       CONSIDERATIONS       AND       EXPECTED       MODEL 
RESPONSE 

Previous  studies  of  isolated  eddies  have  examined  in  detail  the  propagation  and 

decay  of  isolated  eddies  in  quiescent  P-plane  backgrounds  (McWilliams  and  Flierl, 

1979;  Mied  and  Lindemann,  1979).    At  mid-latitudes,  isolated  eddies  have  significant 

propagation  tendencies  associated  with  the  planetary'  vorticity  gradient  (P)  and  can 

decay  through  Rossby  wave  radiation  and  frictional  processes.    These  studies  showed 

that  westward  propagation  speed  is  limited  by  the  Rossby  long  wave  speed: 

cmax  =  -PRd2  (eqn2.U) 

The  value  of  cmax  at  low  latitudes  is  approximately  4  km/day;  however,  at  80 °N, 
cmax  1S  on  ^e  orc*er  °f  -01  km/day.  The  decay  time  associated  with  Rossby  wave 
radiation  (Smith  and  Reid,  1982)  is  substantially  longer  at  high  latitudes.  Simulations 
reported  here  focus  on  the  short  time  interaction  (less  than  1  week)  of  isolated  eddies 
and  ice.  It  is  anticipated  that  eddy  evolution  and  decay  for  these  short  time  cases  are 
related  wind  or  to  ice  processes  and  not  to  the  aforementioned  processes,  which  occur 
on  longer  time  scales. 

The  previous  wind  forced  marginal  ice  zone  ocean  models  give  an  indication  of 
the  expected  Ekman  response  in  the  ocean  and  ice  edge.  Hakkinen  (1986a,  b)  shows 
that  10  m/s  along-ice-edge  winds  give  an  upwelled  or  downwelled  interface  amplitude 
of  ~10  m  and  an  ice  edge  jet  velocity  of  9  cms  after  a  10-15  day  spin-up  period. 
Additionally,  ice  banding  simulations  require  at  least  one  wind  reversal  to  occur.  As 
simulations  here  are  integrated  for  shorter  periods  with  weaker  winds,  weak  ice  edge  jet 
velocities  and  no  ice  banding  are  expected. 
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III.  EXPERIMENTS 

A.       WIND  FORCING  DIRECTIONS 

An  understanding  of  the  effect  of  wind  direction  on  ice  motion  is  the  goal  of  the 
following  experiments.  Variable  in  the  experiments  are  the  wind  magnitude,  direction, 
and  duration  for  a  given  direction.  Values  for  wind  direction  are  measured 
counterclockwise  in  degrees  from  the  positive  x-axis  of  the  model  domain.  Since  the 
model  domain  is  rotated  30  counterclockwise,  a  model  input  wind  direction  of  025° 
means  that  on  the  globe,  the  wind  would  be  generally  from  the  southwest.  This 
correlates  to  a  downwelling  favorable  wind  condition  when  viewed  relative  to  the  initial 
concentrated  ice  band.    Table  1  lists  the  wind  directions  used  in  simulations. 


TABLE  1 

WIND  DIRECTIONS  USED  IN  SIMULATIONS 

RELATIVE 
WIND 

MODEL 
INPUT  WIND 

TRUE 
WIND 

Upwelling  Favorable 

205° 

235° 

Downwelling  Favorable 

025° 

055° 

On-ice 

115° 

145° 

Off- ice 

295° 

325° 

A  model  input  wind  direction  of  205°  would  be  a  true  wind  generally  from  the 
northeast;  this  correlates  to  an  upwelling  favorable  wind.  Winds  perpendicular  to  the 
downwelling  or  upwelling  wind  directions  are  considered  to  be  generally  on-ice  and  off- 
ice winds.  Along-ice-edge  wind  forcing  (0°  and  180°)  induces  a  cross-ice-edge  Ekman 
drift  in  the  ice.  The  025°  and  205°  angles  chosen  by  Smith  et  al.  (1987)  are  intended 
to  minimize  (but  not  eliminate)  the  cross-ice-edge  motion,  thereby  maintaining  the 
concentrated  ice  band  near  mid-domain  and  thus  reducing  boundary  effects.    The  most 
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important  aspect  of  these  experiments  is  not  the  true  direction  of  the  wind,  but  the 
relative  direction  of  the  wind  in  relation  to  the  concentrated  ice  band. 

B.  PRELIMINARY  EXPERIMENTS 

Before  proceeding  to  the  more  realistic  wind-driven  simulations,  the  response  of 
the  ice  edge  to  ocean  eddies,  without  wind  forcing,  is  considered  first. 
1.  Simulations  with  No  Wind 

a.  Cyclonic  Eddy  Initialization 

Figure  3.1  shows  the  evolution  of  the  concentrated  ice  band  under  the 
influence  of  a  cyclonic  ocean  eddy  without  wind  forcing.  By  day  one,  the  eddy  is 
affecting  the  ice  band  in  such  a  way  that  the  ice  band  contains  a  sinuous  pattern  that 
by  day  two  is  more  pronounced.  By  day  three  the  ice  clearly  exhibits  a  circular, 
cyclonically  turning  eddy,  which  is  completely  distinct  from  the  original  ice  band  by 
day  four. 

b.  Anticy clonic  eddy  initialization 

As  would  be  intuitively  expected,  the  anticyclonic  case  (Figure  3.2)  with  no 
wind  forcing  exhibits  nearly  the  same  result  as  above,  except  that  the  sense  of  rotation 
is  opposite  from  the  cyclonic  case. 

The  results  indicate  that  the  ice  velocity  evolves  rapidly  to  the  ocean 
velocity  and  is  simply  advected  radially.  These  results  are  consistent  with  Smith  el  al. 
(1987),  where  a  broader  ice  concentration  band  was  considered.  Smith  et  al.  (1987) 
show  that  the  ice  momentum  equation  in  no  wind  simulations  equilibrates  to  a 
gradient  balance.  For  the  cyclone  case,  a  radially  inward  pressure  force  in  the  ice 
associated  with  the  ocean  surface  slope  is  balanced  by  the  radially  outward  Coriolis 
and  centrifugal  forces.  Likewise,  for  an  anticyclone  case,  the  radially  outward  pressure 
gradient  and  centrifugal  forces  are  balanced  by  the  inward  Coriolis  force.  The  ice 
momentum  balance  was  shown  to  equilibrate  rapidly  in  less  than  one  day. 

C.  BASE  CASE 

Numerous  experiments  were  conducted  using  the  model  previously  discussed.  To 
assess  the  model  output  variations  resulting  from  parameter  changes,  one  simulation  is 
chosen  to  be  a  base  case.  This  simulation  was  selected  as  the  reference  case  because  of 
its  resemblance  to  the  observations  of  Johannessen  et  al.  (1987).  A  summary  of  the 
model  parameters  chosen  for  the  base  case  is  listed  in  Table  2. 
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Figure  3.1     Preliminary  Experiment  (Cyclone). 


Evolution  of  (a)  ice  concentration,  (b)  ocean  upper  laver  potential  vorticitv,  (c)  ocean 
lower  layer  potential  vorticitv,  (dj  ocean  surface  heieht  anomalv,  and"  (e)  ocean 
interface  height  anomaly.   Contour  interval  for  ice  concentration  is  0.t)5. 
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Figure  3.2    Preliminary  Experiment  (Anticyclone). 

Evolution  of  (a)  ice  concentration,  (b)  ocean  upper  layer  potential  vorticity,  (c)  ocean 
lower  laver  potential  vorticity,  (dj  ocean  surlace  height  anomalv.  ana  (e)  ocean 
interface  height  anomaly.    Contour  interval  for  ice  concentration  is  0.1)5. 
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TABLE  2 

BASE  CASE  MODEL  PARAMETERS 

PARAMETER 

SYMBOL 

MODEL  INPUT 

Initial  Upper  Layer  Thickness 

Hl 

50  m 

Initial  Lower  Layer  Thickness 

H2 

4000  m 

Initial  Ice  Thickness 

D 

2  m 

Initial  Ice  Band  Concentration 

A 

.05  to  .75 

Density  of  Ice 

Pi 

9l0kg/m3 

Grid  Spatial  Resolution 

Ax 

1000  m 

Time  Increment 

At 

600  s 

Bottom  Topography 

Flat 

Eddy  Radius 

5  km 

Eddy  Rotation 

Cyclonic 

Wind  Speed 

3  m,  s 

Wind  Direction 

205° 

1.   Experiment  No.  1  (Base  Case:  Cyclone,  Light  Upwelling  Winds) 

Using  the  input  parameters  listed  in  Table  2,  an  entire  range  of  experiments 
can  be  discussed  and  compared  to  this  reference  case.  Light  winds  (3  m/s)  were  chosen 
because  they  are  often  observed,  but  until  now  have  not  been  included  in  model 
simulations. 

With  constant  light  upwelling  favorable  winds  (3  m/s),  the  concentrated  ice 
band  responds  to  the  cyclonic  motion  of  the  ocean  eddy  for  the  first  day  (see  Figure 
3.3).  However,  by  day  two,  the  effect  of  the  wind  causes  a  strong  projection  to  the 
south  of  the  ice  band.  This  projection  of  more  concentrated  ice  becomes  an  enclosed 
cyclonic  feature  by  day  three.  By  day  four  the  cyclonic  feature  has  nearly  broken  away 
from  the  original  concentrated  ice  band.  As  a  result  of  the  growth  of  the  projection 
first  noted  at  day  two,  there  is  a  pocket  of  more  concentrated  ice  that  is  cut  off  from 
ice  of  equal  concentration  that  has  formed  to  the  south.    Also  evident  in  the  figure  is  a 
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significant  northward  Ekman  drift  of  the  ice  edge;  the  resulting  signature  of  the  ocean 
eddy  in  the  ice  is  dependent  upon  cross-ice-edge  motion.  Figure  3.3  also  shows  that  3 
m/s  winds  produce  no  noticeable  effect  on  the  upper  ocean  cyclone.  Figure  3.3  also 
shows  that  insignificant  upwelling  occurs  for  this  wind  magnitude.  The  interface 
anomaly  indicates  downwelling  at  the  top  boundary  of  the  model  domain  and 
upwelling  at  the  bottom  boundary  of  the  model  domain  as  expected  from 
considerations  of  Ekman-induced  motion  at  a  boundary.  Figure  3.4  shows  the  ice 
pattern  concentration  at  day  four  and  day  five. 

Figure  3.5  shows  how  the  velocity  of  the  ice  is  affected  by  the  wind  and  the 
ocean  eddy  for  the  base  case.  The  ice  toward  the  upper  (northern)  boundary,  over  the 
ocean  eddy,  is  moving  faster  than  the  ice  located  toward  the  bottom  (southern) 
boundary,  over  the  ocean  eddy.  The  northern  ice  moves  faster  because  the  ocean  eddy 
and  the  wind  are  working  together;  in  the  southern  part  of  the  model  domain,  the 
ocean-eddy-driven  component  and  the  wind-driven  component  of  ice  motion  are 
opposed. 

The  ice  momentum  balance  for  10  m/s  wind  driven  cases  is  shown  in  Smith 
et  al.  (1987).  The  gradient  ice  momentum  balance  seen  above  in  no  wind  cases 
becomes  a  secondary  balance.  The  dominant  force  balance  in  the  ice  becomes  that 
between  Tai,  the  air-ice  stress  in  the  direction  of  the  wind,  and  TIW,  the  ice-water  dran. 
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Figure  3.3     Experiment  No.  1  (Base  Case). 

Evolution  of  (a)  ice  concentration,  (b)  ocean  upper  layer  potential  vorticity,  (c)  ocean 
lower  laver  potential  vorticity,  (dj  ocean  surface  height  anomaly,  ana"  (e)  ocean 
interface 'height  anomaly.   Contour  interval  for  ice  concentration  is  0.D5. 
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Figure  3.4    Experiment  No.  1  (Base  Case). 

Ice  concentration  at  (a)  dav  4  and  (b)  day  5.   Contour  interval  for  ice  concentration  is 
0.05. 
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Figure  3.5    Plots  of  Ice  Velocity  (Base  Case). 
Ice  velocity  vectors  at  day  3.    Arrows  represent  ice  velocity.    Contours  highlight  areas 


ice  velocity  vectors  at  ai 
of  high  or  low  velocities. 
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D.       VARIATION  FROM  THE  BASE  CASE 

Table  3  lists  the  model  parameters  that  are  varied  in  the  experiments,  beginning 
with  Experiment  No.  1  (base  case).  Since  the  observations  of  Johannessen  et  al.  (1987) 
indicate  that  most  of  the  mesoscale  eddies  observed  in  the  East  Greenland  Current 
MIZ  are  cyclonic,  most  of  the  experiments  are  for  cyclonic  eddies. 


TABLE  3 

MODEL  PARAMETERS  VARIED  IN  EXPERIMENTS 

EXP 

OCEAN 
EDDY 

WIND 
DIRECTION 

WIND 
SPEED 

WIND 
SHIFT 

1 

cyclonic 

205° 

3  m/s 

no 

2 

anticyclonic 

205° 

3  m  s 

no 

3a 

cyclonic 

025° 

3  m/s 

no 

3b 

cyclonic 

025° 

5  m/s 

no 

4 

anticyclonic 

025° 

3  m/s 

no 

5a 

cyclonic 

205°/25° 

10  m/s 

24  hr 

5b 

cyclonic 

205°/25° 

10  m/s 

36  hr 

6a 

cyclonic 

205°/115° 

3  m  s 

24  hr 

6b 

cyclonic 

205°/295° 

3  m/s 

24  hr 

1.  Experiment  No.  2  (Anticyclone,  Light  Upwelling  Winds) 

This  experiment  was  conducted  to  observe  the  effect  on  the  concentrated  ice 
band  of  an  anticyclonic  ocean  eddy  in  an  environment  similar  to  that  of  the  base  case. 
The  anticyclonic  eddy  exhibits  no  development  of  a  distinct  cut-off  feature  (Figure  3.6) 
as  in  the  cyclonic  base  case  with  3  m/s  winds.  The  sinuous  pattern  in  the  concentrated 
ice  band  gives  evidence  of  anticyclonic  motion  at  day  one,  which  becomes  pronounced 
at  dav  two,  but  bv  dav  three  the  anticvclonic  signature  is  decreasing.  Bv  dav  four  the 
ice  pattern,  which  appears  distinctly  different  from  previous  cases,  is  weakly 
anticyclonic  but  by  no  means  representative  of  the  anticyclonic  eddy  that  is  beneath 
the  ice.  Smith  et  al.  (1987)  found  that  10  m/s  wind  in  this  direction  could  destroy  the 
upper  ocean  cyclone. 
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Figure  3.6  shows  that  there  is  no  great  change  in  the  upper  ocean  vorticity 
pattern  or  upper  layer  pressure,  indicating  that  3  m/s  winds  are  of  insufficient  strength 
to  alter  the  upper  ocean  eddy.  However,  this  experiment  shows  that  3  m/s  winds  have 
a  significant  effect  on  the  ice.  This  case  also  shows  the  cross-ice-edge  drift  of  the 
concentrated  ice  band  to  the  north.  Figure  3.7  shows  how  the  velocity  of  the  ice  is 
affected  by  the  wind  and  the  ocean  eddy  for  this  case.  In  this  case  the  wind  and  ocean 
are  opposing  effects,  which  results  in  the  relatively  flat  ice  pattern  by  day  three. 
2.  Experiments  No.  3a,  3b  (Cyclone,  Downwelling  Winds) 

a.  Case  for  3  m/s  Hinds 

Experiment  3a  (Figure  3.8)  was  designed  to  investigate  the  change  in  ice 
pattern  for  winds  from  a  direction  opposite  from  the  base  case.  The  initial  conditions 
for  this  experiment  are  the  same  as  Experiment  No.  1  (base  case),  with  the  exception  of 
the  wind  direction,  which  is  changed  180°  from  the  base  case  so  that  a  downwelling 
favorable  wind  is  used.  With  these  opposite  direction  winds,  no  distinct  cut-off  feature 
develops,  unlike  the  base  case.  At  no  point  in  time  is  there  a  concentration  of  ice  over 
the  eddy  similar  to  that  found  in  the  base  case.  By  day  five  there  is  little  evidence  from 
the  ice  pattern  of  a  cyclonic  ocean  eddy  under  the  ice,  even  though  the  potential 
vorticity  of  the  ocean  eddy  remains  nearly  at  the  same  value  as  that  at  day  zero.  Note 
also  that  the  cross-ice-edge  drift  is  now  to  the  south,  opposite  from  that  for  upwelling 
favorable  wind  conditions.  The  southward  movement  of  the  concentrated  ice  band 
creates  a  situation  that  is  unfavorable  for  any  distinct  ice  pattern  to  break  away  from 
the  ice  band. 

The  differences  in  the  upwelling  and  downwelling  cases  can  be  understood 
from  vorticity  considerations.  The  air-ice  stress  tai  is  weighted  by  ice  concentration,  A, 
in  the  ice  momentum  equations.  Therefore,  moving  from  high  concentration  toward 
low  concentration,  the  air-ice  coupling  becomes  weaker.  Upwelling  winds  thus  provide 
cyclonic  vorticity  to  the  ice  and  downwelling  winds  provide  anticyclonic  vorticity  to  the 
ice.  The  winds  thus  augment  the  ocean  induced  ice  motion  in  the  cyclone  case  and  are 
opposed  to  the  ocean  induced  ice  motion  in  the  anticyclone  case. 

b.  Case  for  5  m/s  Hinds 

Case  3b  (Figure  3.9)  is  the  same  as  the  above  case  with  the  wind  increased 
to  5  m/s  from  the  previous  3  ms.  This  small  increase  shows  how  much  quicker  the 
response  of  the  ice  pattern  is  to  slightly  greater  wind  forcing.  Other  experiments,  not 
shown  here,  show  that  as  the  wind  speed  continues  to  increase  for  the  downwelling 
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favorable  direction,  the  ice  pattern  responds  sooner.  This  illustrates  that  the  wind 
begins  to  dominate  the  ice  momentum  balance  over  the  ocean  forcing  as  winds  increase 
from  3  m/s.  Figure  3.9  shows  that  5  ms  winds  begin  to  induce  a  downwelled  response 
at  the  ice  edge  near  the  end  of  5  days.  The  upper  ocean  pressure  field  indicates  a 
weakening  of  the  ocean  cyclone.  Also  evident  is  an  increased  cross-ice-edge  motion  to 
the  south  of  the  model  domain. 

3.  Experiment  No.  4  (Anticyclone,  Downwelling  Winds) 

This  experiment  (Figure  3.10)  was  done  to  observe  the  effect  of  an  anticyclone 
in  the  presence  of  downwelling  favorable  winds.  Based  on  the  results  of  Smith  et  al. 
(1987),  the  expections  were  that  the  anticyclonic  signature  in  the  ice  would  not  be 
destroyed  as  quickly  as  in  Experiment  No.  2  where  upwelling  favorable  winds  were 
modeled.  The  anticyclonic  signature  is  more  pronounced  because  the  vorticity  input  to 
the  ice  from  the  ocean  and  air  are  both  anticyclonic;  in  Experiment  No.  2  the  vorticity 
input  to  the  ice  form  the  ocean  and  air  are  opposed. 

For  the  first  two  days  of  this  simulation,  the  anticyclonic  signature  in  the  ice 
develops  much  the  same  as  in  the  cyclonic  simulation  (Experiment  3a),  but  with  an 
opposite  sense  of  rotation.  The  difference  between  this  simulation  and  the  base  case  is 
that  instead  of  having  a  higher  ice  concentration  over  the  ocean  eddy  center,  a  region 
of  less  concentrated  ice  becomes  centered  over  the  ocean  eddy.  Another  difference  is 
the  shape  of  the  lower  ice  edge  of  the  concentrated  ice  band.  In  this  case  (Figure  3.1 1) 
the  lower  ice  edge  is  smoother  and  less  remarkable  than  found  in  the  base  case. 

By  the  above  arguments  of  vorticity  input  to  the  ice  induced  by  the  wind,  this 
case  should  be  the  mirror  image  of  the  base  case.  A  comparison  of  Figure  3.10  to 
Figure  3.3  shows  this  not  to  be  the  case.  The  explanation  for  his  lack  of  symmetry  lies 
in  the  Ekman  induced  cross-ice-edge  component  of  motion.  In  the  base  case,  the  ice 
band  drifted  northward,  away  from  the  stationary  ocean  eddy.  In  this  simulation  the 
ice  band  has  a  southward  drift,  bringing  higher  ice  concentration  across  the  ocean 
eddy. 

4.  Experiments  No.  5a,  5b  (Cyclone,  180°  Reversal  in  Winds) 

For  comparison  with  the  results  of  Smith  et  al.  (1987),  where  10  ms  constant 
direction  winds  are  considered,  wind  reversals  are  included  next. 

a.    Case  for  10  m/s  Hinds  with  24  Hour  Reversal 

This  experiment  (Figure  3.12)  shows  the  results  of  10  m/s  winds  being 
instantaneously  reversed   from  upwelling   favorable  winds   to   downwelling   favorable 
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winds  every  24  hours.  This  was  repeated  for  a  total  of  six  days.  Using  this  technique, 
the  ice  showed  much  the  same  response  as  the  case  with  no  wind;  however,  the  ice  over 
the  eddy  is  somewhat  more  loosely  packed.  This  experiment  shows  the  vast  difference 
in  the  ice  pattern  by  including  a  wind  shift  in  the  model.  By  applying  the  results  of 
Experiment  No.  3b,  a  10  ms  downwelling  wind  would  quickly  dominate  the  ice  pattern 
initially  induced  by  the  ocean  eddy.  However,  by  daily  reversals  of  the  direction  of  the 
wind,  the  wind  works  in  conjunction  with  the  eddy  to  form  a  concentrated  ice  pattern 
over  the  ocean  eddy.  In  terms  of  vorticity,  the  wind  is  alternately  providing  cyclonic 
(during  upwelling)  and  anticyclonic  (during  downwelling)  relative  vorticity  to  the  ice. 

This  experiment  was  also  repeated  using  a  number  of  other  wind  speeds. 
The  results  were  similar;  for  each  case  with  wind  speeds  less  than  10  ms,  ice  spirals 
radially  over  the  eddy  when  the  wind  is  shifted  every  24  hours.  This  is  in  contrast  to 
the  results  of  Smith  et  al.  (1987),  where  a  10  m,  s  wind  was  shown  to  straighten  the  ice 
edge,  overcoming  the  effect  of  the  ocean  eddy.  In  Smith  et  al.  (1987)  it  was  also  shown 
that  constant  direction  upwelling  (downwelling)  winds  destroyed  an  upper  ocean 
cyclone  (anticyclone).  With  wind  reversals  of  24  hours,  this  eddy  destruction 
mechanism  no  longer  is  observed.  This  is  consistent  with  a  vorticity  argument 
provided  by  Smith  et  al.  (1987).  As  air-ice  coupling  is  stronger  than  air-water  coupling 
by  virtue  o[  the  drag  coefficients,  the  momentum  input  to  the  ocean  by  winds  is  larger 
under  the  ice  than  in  the  open  ocean.  Upwelling  (downwelling)  winds  hence  provide 
cyclonic  (anticyclonic)  vorticity  to  the  upper  ocean.  In  24  hour  wind  reversals  here,  the 
wind  alternately  provides  cyclonic  and  anticyclonic  relative  vorticity,  thus  allowing  the 
ocean  eddy  to  persist.  Figure  3.13  shows  the  vorticity  in  the  ice  for  a  concentrated  ice 
band  with  no  eddy  in  the  presence  of  upwelling  and  downwelling  winds. 
b.    Case  for  10  mjs  Winds  with  36  Hour  Reversal 

All  other  parameters  are  the  same  as  the  previous  case,  except  that  the 
wind  shifts  at  36  hour  intervals.  The  difference  in  the  ice  pattern  evolution  is  very 
striking.  This  case  (Figure  3.14)  does  not  develop  a  concentrated  ice  formation  over 
the  ocean  eddy.  The  concentrated  ice  band  shape  is  irregular  throughout  the 
simulation.  The  ice  pattern  gives  no  strong  evidence  of  an  ocean  eddy,  even  though 
the  ocean  eddy  persists  during  this  period. 
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5.  Experiments  No.  6a,  6b  (Daily  90°  Shift  in  Light  Winds) 

a.  Shift  from  Upwelling  to  On-Ice  Winds  Kith  a  Cyclonic  Eddy 

Here  the  wind  was  chosen  as  in  the  base  case  with  a  wind  speed  of  3  m's. 
In  this  experiment  (Figure  3.15)  the  wind  was  shifted  90°  back  and  forth  every  24 
hours,  so  that  the  wind  would  shift  from  upwelling  favorable  to  on-ice  after  24  hours, 
then  shift  back  after  the  next  24  hours.  The  wind  shift  was  continued  for  four  days. 
The  ice  pattern  develops  a  strong  cyclonic  cut-off  feature  at  the  eastern  edge  of  the 
concentrated  ice  band  that  is  much  more  pronounced  than  the  cut-off  feature  to  the 
bottom  of  the  model  domain  found  in  any  previous  experiment.  The  wind  from  the 
southeast  is  an  on-ice  wind,  which  drives  the  concentrated  ice  band  to  the  top  of  the 
model  domain.  The  effect  on  the  ice  of  upwelling  winds  with  a  cyclonic  ocean  eddy 
alone  creates  a  near-break-away.  The  addition  of  on-ice  winds  greatly  increases  the 
concentration  of  the  southward  projection  away  from  the  concentrated  ice  band 
(Figure  3.16).  The  enhancement  of  cross-ice-edge  motion  by  shifting  the  wind  shows 
how  sensitive  the  response  of  the  ice  is  to  the  wind  direction. 

b.  Shift  from  Upwelling  to  Off-Ice  Hinds  with  a  Cyclonic  Eddy 

In  this  case  the  wind  is  shifted  from  upwelling  favorable  to  off-ice  winds 
every  24  hours.  Figure  3.17  shows  that  the  amount  of  ice  over  the  ocean  eddy  in  this 
case  is  less  than  in  the  on-ice  case.  This  is  due  to  the  competing  forces  driven  by  the 
wind  and  the  ocean  eddy. 

A  comparison  to  the  results  of  this  simulation  with  the  base  case  illustrates 
an  important  point.  The  breaking  away  of  ice  from  the  ice  edge  is  dependent  upon  the 
ice  edge  retreating  away  from  the  eddy.  In  this  simulation,  no  cross-ice-edge  motion 
occurs  and.  hence,  no  break  away  of  ice  from  the  edge  occurs.  Thus,  the  process  of  ice 
trapping  over  eddies  is  sensitive  to  wind  direction  and  the  resulting  cross  ice  edge  shifts. 
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Figure  3.6    Experiment  No.  2  (Anticyclone,  Upwelling  Winds). 

Evolution  of  (a)  ice  concentration,  (b)  ocean  upper  laver  potential  vorticity,  (c)  ocean 
lower  laver  potential  vorticitv,  (d)  ocean  surface  height  anomaly,  ancf  (e)  ocean 
interface 'height  anomalv.   Contour  interval  for  ice  concentration  is  0.D5. 
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Figure  3.7    Experiment  No.  2  (Anticyclone,  Upwelling  Winds). 

Ice  Velocity  Vectors  at  day  3.   Arrows  represent  ice  velocity.   Contours  highlight  areas 
of  high  or  low  velocities. 
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Figure  3.8     Experiment  No.  3a  (Cyclone,  Light  Downwelling  Winds). 

Evolution  of  (a)  ice  concentration,  (b)  ocean  upper  layer  potential  vorticity,  (c)  ocean 
lower  laver  potential  vorticity,  (d)  ocean  surface  height  anomaly,  and  (e)  ocean 
interface  height  anomalv.   Contour  interval  for  ice  concentration  is  0.05. 
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Figure  3.9    Experiment  No.  3b  (Cyclone,  5  m/s  Downwelling  Winds). 

Evolution  of  (a)  ice  concentration,  (b)  ocean  upper  laver  potential  vorticitv,  (c)  ocean 
lower  layer  potential  vorticity,  (dj  ocean  surlace  height  anomalv,  ancf  (e)  ocean 
interface  height  anomaly.   Contour  interval  for  ice  concentration  is  0.05. 
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Figure  3.10    Experiment  No.  4  (Anticyclone,  Light  Downwelling  Winds). 


Evolution  of  (a)  ice  concentration,  (b)  ocean  upper  laver  potential  vorticity,  (c)  ocean 
lower  laver  potential  vorticitv,  (d)  ocean  suriace  height  anomaly,  ana  (e)  ocean 
interface 'height  anomaly.   Contour  interval  for  ice  concentration  is  0.05. 
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Figure  3.11     Experiment  No.  4  (Anticyclone,  Downwelling  Winds). 

Ice  concentration  at  (a)  day  4  and  (b)  day  4.5.    Contour  interval  for  ice  concentration 
is  0.05. 
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Figure  3.12    Experiment  No.  5a  (Cyclone,  Daily  180°  Reversal  in  10  m/s  Winds  ). 

Ice  concentration  at  (a)  day  3  (b)  day  4  and  (c)  day  5.    Contour  interval  for  ice 
concentration  is  0.05. 
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Figure  3.13     Vorticity  Input  to  the  Ice. 

Vorticity  across  the  ice  band  for  (a)  downwelling  winds  and  (b)  upwelling  winds 
Where  J  is  the  y-coordinate  grid  point. 
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Figure  3.14    Experiment  No.  5b  (Cyclone,  36  hr  180°  Reversal  in  10  m/s  Winds). 

Ice  concentration  at  (a)  day  1  (b)  day  2  and  (c)  day  3.    Contour  interval  for  ice 
concentration  is  0.05. 
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Figure  3.15    Experiment  No.  6a  (Cyclone,  Daily  90°  Shift  to  Light  On-Ice  Winds). 

Evolution  of  (a)  ice  concentration,  (b)  ocean  upper  laver  potential  vorticitv,  (c)  ocean 
lower  laver  potential  vorticity,  (d)  ocean  surface  height  anomalv,  and  (e)  ocean 
interface 'height  anomaly.   Contour  interval  for  ice  concentration  is  0.D5. 
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Figure  3.16    Experiment  No.  6a  (Cyclone,  Daily  90°  Shift  to  Light  On-Ice  Winds). 
Ice  concentration  at  day  4.   Contour  interval  for  ice  concentration  is  0.05. 
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Figure  3.17    Experiment  No.  6b  (Cyclone,  Daily  90°  Shift  to  Light  Off-Ice  Winds). 

Evolution  of  (a)  ice  concentration,  (b)  ocean  upper  layer  potential  vorticity,  (c)  ocean 
lower  laver  potential  vorticitv,  (d)  ocean  suriace  height  anomalv,  and  (e)  ocean 
interface 'height  anomaly.   Contour  interval  for  ice  concentration  is  0.05. 
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IV.  DISCUSSION  AND  CONCLUSIONS 

A.  THE  IMPORTANCE  OF  CYCLONIC  EDDIES  IN  THE  MIZ 

Johannessen  et  al.  (1987)  found  that  most  of  the  mesoscale  eddies  in  the 
Marginal  Ice  Zone  of  the  Cast  Greenland  Current  were  cyclonic  in  nature.  Only  2  of 
the  15  eddies  observed  in  Johannessen  et  al.  (19S7)  were  anticyclones.  The  authors 
hint  that  the  two  anticyclones  may  be  the  result  of  topography  because  these 
anticyclonic  eddies  were  found  only  over  topographic  highs.  The  results  of  these 
simulations  find  that  anticyclonic  eddies  interacting  with  an  ice  edge  do  not  leave  as 
clear  of  a  signature  of  their  presence  at  the  ice  edge.  Simulated  cyclonic  eddies 
interacting  with  the  ice  edge  leave  their  signature  for  much  longer  periods  and,  in  some 
special  circumstances,  enhance  the  signature  of  the  ocean  eddy  by  aiding  in 
concentrating  the  ice  pattern  over  the  eddy. 

B.  APPLICABILITY  OF  WIND  SHIFT 

The  ice  banding  discussed  by  Hakkinen  (1986a)  was  not  observed  as  expected  for 
two  reasons: 

1.  The  simulations  in  this  study  were  not  run  out  to  the  time  for  which  ice 
banding  was  observed  in  Hakkinen  (1986a). 

2.  The  effects  of  ice  banding  could  be  masked  by  the  interference  caused  by  the 
ocean  eddy. 

Hakkinen  (1986a)  chose  to  vary  winds  sinusoidally  with  a  period  of  four  days  to 

resemble  the  passage  of  atmospheric  cyclones  over  the  ice  edge.    This  is  an  adequate 

approximation  for  the  scale  that  Hakkinen  (1986a)  was  simulating. 

Shifting  the  winds  180°  at  intervals  of  24  hours  seemed  to  have  little  effect  on 

the  modeled  ice  pattern  after  a  few  days.    However,  the  cases  where  the  wind  was 

shifted  90°  developed  very  interesting  ice  patterns.    The  striking  resemblance  of  these 

modeled  ice  patterns  to  those  that  have  been  observed  in  periods  of  moderate  winds 

leads  to  the  conclusion  that  the  wind  is  extremely  important  in  the  development  of 

these  ice  patterns  in  the  MIZ. 

C.  NEAR-BREAKAWAY  ICE  FORMATION 

Another  interesting  result  from  the  numerical  simulations  is  that  no  simulated 
eddy  ever  broke  completely  away  from  the  concentrated  ice  edge.    There  were  some 
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cases  where  the  concentration  in  the  center  of  a  cyclonic  ice  pattern  was  considerably 
greater  than  the  ice  concentration  surrounding  it,  but  not  quite  a  complete  separation. 
Thus,  it  can  be  concluded  that  the  very  interesting  eddies  observed  during  previous 
studies  can  not  be  explained  merely  by  the  existance  of  an  eddy  at  the  ice  edge  and  by 
a  regular  pattern  of  wind  shifts  as  used  in  this  study.  This  suggests  that  the  physical 
processes  leading  to  a  cut-off  ice  eddy  are  more  complex  than  can  be  simulated  in  this 
study. 

D.       CONCLUSIONS 

Observations  indicate  the  presence  of  large  (60  km)  and  mesoscale  (5  to  10  km) 
eddies  in  the  East  Greenland  Current  Marginal  Ice  Zone.  This  study  investigates  the 
mesoscale  eddy  interaction  with  the  ice  edge.  Previous  model  studies  of  ocean  eddy 
interaction  with  the  ice  edge  considered  only  strong  (10  m/s)  and  constant  direction 
winds.  This  study  is  concerned  with  the  periods  where  the  wind  speed  is  considered 
light  and  variable.  The  light  winds  and  24  hour  interval  wind  shift  are  viable 
considerations  because,  as  shown  in  Figure  1.2  and  Figure  1.3,  there  are  significant 
periods  where  the  winds  are  light,  and  periods  where  the  wind  direction  shifts  several 
times  over  short  periods.  These  short  duration  wind  shift  periods  may  be  because  of 
possible  local  atmospheric  effects  that  are  of  a  scale  smaller  than  the  circulation 
pattern  considered  by  Hakkinen  (1986a). 

Smith  et  al.  (1987)  have  shown  the  effect  of  eddies  in  the  presence  of  no  wind 
and  in  the  case  of  strong  constant  direction  (10  m/s)  winds  and  a  broad  ice  band.  The 
experiments  completed  here  have  verified  the  results  of  the  no  wind  case  for  the 
cyclonic  and  anticyclonic  eddies  simulated  by  Smith  et  al.  (1987)  for  a  narrow  marginal 
ice  zone.    Specifically,  the  experiments  show  that: 

1.  Three  m/s  winds  are  too  weak  to  induce  upwelling  or  downwelling  at  the  ice 
edge  in  less  than  of  6  days. 

2.  The  effect  of  3  m/s  winds  on  ice  motion  are  comparable  to  the  effect  of  10  cm,  s 
ocean  eddy  motion. 

3.  Upwelling  winds  provide  cyclonic  vorticity  to  the  ice  field  and,  thus,  augment 
cyclonic  ocean-eddy-induced  ice  motion.  Conversely,  downwelling  winds 
provide  anticyclonic  vorticity  to  the  ice,  countering  cyclonic  ocean-eddy-induced 
ice  motion. 

4.  Cross-ice-edge  motion  is  very  sensitive  to  magnitude  and  direction  of  the  wind, 
and  the  resulting  signature  of  the  ocean  eddy  in  the  ice  is  sensitive  to  cross-ice- 
edge  motion. 
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5.  Ice  can  become  trapped  over  an  ocean  eddy,  and  can  nearly  break  away  from 
the  ice  edge  as  winds  cause  cross-ice-edge  motion  of  the  MIZ  to  the  north. 

6.  Alternating  upwelling  downwclling  10  m  s  winds  provide  alternating 
cyclonic/anticyclonic  relative  vorticity  to  the  ice.  If  the  period  of  reversal  is 
short  (24  hours)  the  spin-up; spin-down  effects  cancel  each  other,  leaving  an  ice 
edge  signature  comparable  to  no  wind  cases.  If  the  duration  is  increased  to  36 
hours,  the  wind  vorticity  input  dominates  over  the  ocean  vorticity  input  and  the 
ice  shows  no  eddy  motion. 

7.  The  near-break-away  of  ice  from  an  ice  edge  occurs  in  these  simulations  only 
when  significant  wind  induced  cross-ice-edge  motion  causes  the  initial  ice  edge 
to  retreat  to  the  north,  away  from  the  ocean  eddy  leaving  trapped  ice  over  the 
eddy.  This  cross-ice-edge  motion  is  sensitive  to  magnitude,  direction,  and 
duration  of  wind  events. 

E.       RECOMMENDED  FOLLOW-ON  STUDIES 

Because  the  region  of  the  East  Greenland  Current  has  been  observed  to  be  the 
major  source  of  heat  exchange  between  the  Atlantic  water  and  arctic  water,  further 
studies  should  investigate  the  addition  of  thermodynamic  effects  on  mesoscale  eddies. 
Additionally,  as  anticyclones  would  melt  ice,  the  response  of  ice  to  anticyclones  and 
cyclones  should  show  interesting  differences. 
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APPENDIX 
SYMBOLS  AND  NOTATION 
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Pt 
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ciw 
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Pl 

P2 
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Ice  concentration 

Eddy  maximum  amplitude 

Laplacian  lateral  friction  coefficient  for  ice  concentration 

Laplacian  lateral  friction  coefficient 

Laplacian  lateral  friction  coefficient  for  ice  mass 

Variation  of  Coriolis  parameter  with  latitude 

Topographic  P 

Air-ice  interfacial  stress  coefficient 

Air-water  interfacial  stress  coefficient 

Ice-water  interfacial  stress  coefficient 

Phase  speed  of  perturbation  wave 

Ice  thickness  distribution 

Grid  spatial  resolution 

Time  increment 

Kronecker  delta  function 

Coriolis  parameter  for  mean  latitude 

Gravitational  acceleration 

Reduced  gravitational  acceleration 

Nondimensional  eddy  size 

Upper  layer  mean  thickness 

Lower  layer  mean  thickness 

Instantaneous  layer  thickness 

e-folding  scale  of  the  ocean  eddy 

Ice  mass 

Pressure  in  the  upper  layer 

Pressure  in  the  lower  layer 

Nondimensional  eddy  strength 

upper,  lower  layer  potential  vorticity 

First  internal  Rossby  radius  of  deformation 

Rossby  number  for  upper(i=  1),  lower(i  =  2)  layers 

Densitv  of  air 


3.8  x  10-l2m"l/s-1 


=  3.2x  lCf9 


-3 


=  2.5x  10 
=  1.4x  10" 

=  7.5x  1CT 


m/p-A 

1.0  km 

600  s 

0  when  i  =  2 

1.43  x  10"4s-1 

9.8  m/s2 

g(P2-Pi)Pi 
L/Rd 

50  m 

4000m 

50  m 


prg'hi 

vmax;PRd 

(f  +  yn, 

[g'H^^H^H^]1/2 

v        TL 

max1J 


1  kg  m 


-3 
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p.  Density  of  ice  =910  kg  m"3 

tai  Air-ice  interfacial  stress  vector 

Taw  Air-water  interfacial  stress  vector 

T1W  Ice-water  interfacial  stress  vector 

t_„_w  duration  of  simulation 

1 1  Id  A 

uw,\'w  Velocities  in  the  x  and  y  directions 

U:,V,  Transport  in  the  x  and  y  directions 

vmax  Open  eddy  maximum  tangential  velocity  =  14  cm  s 

x,y  Cartesian  coordinates  directed  NE  and  NW  respectively 

C  Upper,  lower  layer  relative  vorticity  =  V   x  v- 

V  Gradient  operator  =d!dx  +  d;dy 

V2  Laplacian  operator  =  d2/dx    +  d2,dy2 
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